February 4, 2008 1:24 



WSPC - Proceedings Trim Size: llin x 8.5in beskin 



1 

ON THE CYLINDRICAL GRAD-SHAFRANOV EQUATION 



V. S. BESKIN* and E. E. NOKHRINA 



00 

o 

O 
(N 

X) 



Oh 
I 

O 



> 



O 
oo 
O 



I.E.Tamm Theoretical Department, P.N.Lebedev Physical Institute, Moscow, Russia 

* E-mail: beskin&lpi.ru 
www. tamm. lpi.ru 



The goal of this presentation is in paying attention to the ID cyUndrical version of the Grad-Shafranov (GS) equation. 
In our opinion, this approach is more rich than classical self-similar ones, and more suitable for astrophysical jets we 
observe. In particular, it allows us describing the central (and, hence, the most energetic) part of the flow. 
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1. Introduction 

An activity of many compact objects - Active Galac- 
tic Nuclei (AGNs), Young Stellar Objects (YSOs), 
microquasars - is associated with the highly coUi- 
mated jets. These jets are thought to be a natu- 
ral outlet of an excess angular momentum of a cen- 
tral object and accreting matter [1]. The latest ob- 
servations indicating the jet rotation in AGNs [2] 
and YSOs [3] support this idea. The most attractive 
model for such outflows is the MHD one [1, 4, 5]. 

Of course, the main question within this model 
is the coUimation itself [4-8]. We assume here that 
collimation is due to a finite external gas and/or 
magnetic pressure [9-11]. Indeed, proposing that it 
is the external magnetic field Bcxt ^ 10^^ G that 
plays the main role in the collimation, we obtain 

1/2 

Tjet ~ -Rin (-Bin/-Bext) • Here r is the distance from 
the rotational axis, and the subscripts 'in' correspond 
to the values in the vicinity of the central object. The 
similar evaluation can be obtained for external pres- 
sure pext ~ B^^JStt. E.g., for YSOs (Bin - lO^G, 
Rin ^o) we obtain rjot 10^^ cm, in agree- 
ment with observational data. Accordingly, for AGNs 
{Bin ~ 10^ G, Rin ~ 10" cm) we have rjet ^ 1 pc. 
It means that the external media may indeed play 
important role in the collimation process. 

The internal structure of cylindrical jets was con- 
sidered both for non-relativistic [12, 13] and rclativis- 
tic [9, 14-19] flows. In particular, it was shown that 
for constant angular velocity of plasma Sip it is im- 
possible to obtain reasonable solution with total zero 
electric current [9] , but it can be constructed if the 
angular velocity vanishes at the jet boundary and if 
the external pressure is not equal to zero [11, 17]. 

Another result was obtained for relativistic and 



non-relativistic cylindrical flows [13, 15, 16, 20] is 
that the poloidal magnetic field Bp has a jet-like form 

' / core 

where rcoro = i'in7in/^^- But this relation corresponds 
to a very slow (logarithmic) growth of the magnetic 
flux function: ^'(r) oc Inr. It means that if the jet 
core contains only a small part of the total magnetic 
flux 5*0, the jet boundary is to locate exponentially 
far from the axis, magnetic field being too weak to 
be in equilibrium with the external media. In what 
follows we'll try to resolve this contradiction. 

Thus, we consider the following model: the flow 
crosses all the critical surfaces while the effects of 
the external media are negligible. It allows us to use 
standard values of integrals of motion. As the super- 
sonic wind expands, its pressure becomes compara- 
ble with the external gas and/or magnetic pressure. 
The interaction of a flow with external media results 
in well collimated jet which can be described by ID 
cylindrical equations. 

2. Basic equations 
2.1. Relativistic flow 

For cylindrical flow one can write down electric and 
magnetic fields as well as the four-velocity of a 
plasma u in standard form 

21 



rJ = 

27rr 



— e^, E: 

rc 



u= -^B + 7(17Fr/c)ey. (3) 
n 

Here n is the concentration in the comoving refer- 
ence frame, and 7^ = u^ + 1 is the Lorentz-factor. In 
other words, it is convenient to express all the values 
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in terms of magnetic flux ^ and total electric current 

/, the angular velocity of plasma and the particle 
to magnetic flux ratio rj being constant on the mag- 
netic surfaces: Op = i^F(^), V = Accordingly, 
the trans-fleld GS equation can be rewritten as [21] 

2 



1 d 

r dr 



r dr 



+ 



\ dr 
327r^ d 



dl^F 

d¥ 

G 



0. 



(4) 



Here G = r'^{E - ^pLf + M'^L'^c^ - M'^r'^E^, the 
Alfvenic factor is ^ = 1 - Qpr'^/c^ - M"^, 



n 



(5) 



is the Alfvenic Mach number, fi = rUpC^ + mp^u is the 
relativistic enthalpy, and the derivative d/d'if acts on 
the integrals of motion only. Finally, the relativistic 
Bernoulli equation Up = — — 1 has a form 



where 

K = r^{e'f{A 



d*y 
d^y 



K 



^2 ^2^4 



2 2 



(6) 



M^)+M^r^E^ -M^L^c", (7) 



and e' = E — QpL. Both equations contain relativis- 
tic integrals of motion 

QfI / 



E{'^) = ^i^irjc^ 



27r 



which, as all other invariants, are to be determined 
from boundary and critical conditions. E.g., for in- 
ner part of a flow ^ <^ with zero temperature 
one can choose Of(^') = Oq, r]{^) = r]o, and 

= winc^ + i(g>) = i^f. (9) 

Multiplying now equation (4) on 2Ad^'/dr and 
using equation (6), one can find [17] 



(e') 



/\2 



2».2 



1 + 



6 r2 



Ar^^'^rj^c^ 



+ 



dAf_ 

dr 

{e'r 



d^- de' 



^2^2^,4 dr d^E" 



Aijprj'^c^ 
2 d* dOp 



dr d* 



„ , , d* 1 dr) 
+ 2^44 -4- 



'dP 
S7 



T 



J dr Tj d^ 
d* ds 
jj, dr d* 



where Cg <C c is the sound velocity, and the entropy 
s = .s('J') is the fifth integral of motion. Together 
with Bernoulli equation (6) it forms the system of 
two ordinary differential equations for Mach number 
aj^d magnetic fiux describing cylindrical rela- 
tivistic jet. Clear boundary conditions are 



*(0) = 0, 

-P(''jet) = -fexti 



(11) 

(12) 



where P = B'^/Sn + p is the total pressure. Deter- 
mining the functions A4^(r) and ^'(r), one can find 
the jet radius rjet as well as the profile of the cur- 
rent J, particle energy, and toroidal component of 
the four-velocity using standard expressions 



/ _ L-VlYr^E/c^ 
2^ ~ l-nlr^/c^ - M^' 

1 {E -VIfL) - M^E 
^ ~ 1 - r!2r2/c2 _ X2 ' 

1 {E - nFL)nFr^/c^ - LM^ 



jirirc 



l-f22r2/c2-A^2 



(13) 
(14) 
. (15) 



2.2. Nonrelativistic flow 

In the nonrelativistic limit electric and magnetic 
fields arc determined by general expressions (2). On 
the other hand, equation (3) can be rewritten as 



V = — B + flFretn, 

Pm 



(16) 



where pm = rupn is the mass density and 7?n(*) is 

nonrelativistic particle to magnetic fiux ratio. Ac- 
cordingly, nonrelativistic fluxes of energy E^ and z 
component of the angular momentum L„ are 



27rC77n 



Vi^r sin 9. 



(18) 



Further, algebraic relations (13)-(15) can be rewrit- 
ten as 



where now 



2^ ~ ^''^ 1-M2 ' 
_ 1 riFr2 - L^M^' 
- - i_^2 



(19) 
(20) 

(21) 
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As a result, nonrelativistic Bernoulli equation 



647r'*77^ \dr J 



(1-^2)2 - 1-X2 ' (22) 

together with nonrelativistic limit of equation (10) 



2e„ -2w + O^r^ - (1 - 



dr 
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2„ 



,d* dcn 
dr d* 
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d-^ df^F 
dr d* 

d* d77n 

?7ii dr d^P 



1 /ap 

9s 
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T 



d* ds 



dr d* 



:^3) 



where Cn = ^ ^^F-^^nj determine the structure of 
nonrelativistic cylindrical flow. 

3. Advantages 

Certainly, the approach under consideration is ID as 
well. For this reason, it has some properties similar 
to another self-similar ones. In particular, one can 
easily check that the singularity on the fast magncti- 
sonic surface is absent. On the other hand, singular- 
ity appears on the cusp surface where the factors in 
front of dA<^/dr in (10) and (23) vanish. Neverthe- 
less, in our opinion, this one-dimensional approach 
has clear advantages in comparison with the stan- 
dard self-similar ones [4-7]. 

First of all, it allows us to use any form of the five 
integrals of motion. Indeed, the self-similarity of a 
flow demands definite dependence of invariants which 
may be not correspond to the real boundary condi- 
tions. E.g., for relativistic self-similar flow the angu- 
lar velocity Op is to have the form Op oc r"^ [22]. 
It does not correspond neither to the homogeneous 
stellar rotation, nor to the Keplerian disk rotation. 
Moreover, this dependence has the singularity at the 
rotational axis. Thus, the standard self-similar ap- 
proach cannot describe the central (and, hence, the 
most energetic) part of the flow. 

Further, classical self-similar approach cannot 
describe the region of electric current closure. Finally, 
for relativistic magnetically dominated flow it is more 
convenient to use flrst-order equation (10) instead of 



second order GS equation for which it is necessary to 
be careful in taking into account small but important 
terms ~ 7~2. Indeed, the force balance equation (10) 
does not contain the leading terms and j x B/c as 
they are analytically removed using Bernoulli equa- 
tion. As a result, as 

|peE+jxB/c| l_ 

|jxB/c| '^T^' ^'^^ 

all the terms in equation (10) are of the same order. 

In particular, in the limit r ^ Vcoie, A^^ ^ 1 
equation (10) can be rewritten in the simple form [11] 



dr \ 



0. 

(25) 

Without the last term oc L^{'i) equation (25) results 
in the conservation of the value H 
r2F?7r2 



H 



— const 



(26) 



was found in [26] for conical magnetic field. It is the 
conservation of H that results in the jet-like solution 
(1). Indeed, as r]{^) w const and Of(^') ~ const in 
the center of a jet, we obtain Ai"^ oc r^. Using now 
the definitions A^^ — ^■j^rj'^^/n and nUp = r]Bp (and 
the condition Up « const fulfilled in the very center 
of a flow), we return to (1). But, as we will see, the 
term containing (which appears to be missed pre- 
viously) can be important [17]. It is this term that 
can change the jet-like structure in relativistic case. 

4. Internal structure of cylindrical jets 
4.1. Relativistic flow 
4.1.1. General properties 

The solution of equations (6) and (10) depends es- 
sentially on the Mach number on the rotational axis 
Ml = M^{0) [11]. For Ml » Ml, where Ml, = 7?„ 

(27) 



M^ = MU1 + 



the poloidal magnetic field corresponding to jet-like 
solution (1). On the other hand, for A^g <€. Ml, 



M^ = Ml[l 



2Ml(T \ B.lJ 



(28) 

Here i?L = c/fiF(0), and cr = il^^o/STr^c^/iTyo is the 
Michel magnetization parameter [23] (7 ~ ct for par- 
ticle dominated flow Wp^it ~ Wem)- It means that 
Bp Rs const, i.e., the solution has no jet-like form. 
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As was already stressed, the solution (27) can- 
not be realized in the presence of the external media. 
Hence, for any finite pressure Pgxt rnagnetic field in 
the center of cylindrical jet Bo = 47r7yoM7in/A^o can- 
not be much smaller than Bjnin = 47r77oM7iii/A^cr- It 
gives 



B„ 



1 



(29) 



where B{Ri^) = ^o/t^RI- 



4.1.2. Central core 

Thus, for external magnetic field i?cxt > -Bmin the in- 
ternal structure of a relativistic jet is to be described 
by relations (28). On the other hand, for Bext < -Bmin 
the core with Bp w Bmin is to be formed in the cen- 
ter of a flow (i.e., for r < 7in^L)- In particular, for 
a-'^B{Ri,) < Bext < Bmin (and for r » 7ini?L) the 
solution can be presented as [18] 



M'^ocr", ^ocr'^. 



(30) 



the sum being a + /3 = 3. E.g., for B^xt = Bmin we 
have a = 1, /3 = 2 (cf. 28), and for Bext = a-^B{Ri^) 
we have a = 2, /3 = 1. 

The results presented above were reproduced 
recently both analytically and numerically. In [18] 
it was shown that ID approximation becomes true 
for paraboloidal outflow at large distances from the 
equatorial plane z ^ a'^/^Ri^ where the flow becomes 
actually cylindrical. Up to the distance z = a^iaRh 
one can use the relations (28), so that the poloidal 
magnetic field does not depend on r. The region 
z > cr7i„i?L corresponds to core- like solution (30). 
Nevertheless, the transverse dimension of a jet re- 
mains parabolic: rjet oc z^/^. Numerically the scal- 
ings (30) were confirmed in [24] . 

Remember that the existence of cylindrical core 
with Tcore ^ JinRh was predicted in many papers [14, 
26], but magnetic flux ^core = Trr^o^gBmin inside the 
core was unknown up to now. As we see, in relativis- 
tic case the central core contains only a small part of 
the magnetic flux: 



core _ , Ti 
(J 



in 



*0 



(31) 



Nevertheless, as /? > 0, such core-like flow can exist 
in the presence of external media. 



4.1.3. Bulk acceleration 

As on the fast magnetosonic surface the bulk plasma 
Lorentz-factor jirp) ~ a^l"^ (and, hence, here 
W^part/Wem ~ o""^/^ < 1) [23, 30], the additional 
particle acceleration is possible for r > rp- Using 
equation (14) and relation a -|- /3 = 3 one can find 
that in all region Bext > a~^B(i?L) {z < a'^Rh for 
parabolic flow) the Lorentz-factor for r > jinRh can 
be determined as 



7 w r/Ri,. 
Accordingly, one can write down [17] 



part 



We, 



B(Bl) 



1/2 



-Bext 

-^B{R 



(32) 



(33) 



It means that for Bext ~ cr~'^B[R\,) [z ~ a R\, for 
parabolic flow) where the transverse jet dimension 
rjet ~ c-Rl almost the full energy transformation 
from the Poynting to particle energy flux can be re- 
alized. In particular, for the particle moving along 
parabolic magnetic fleld line one can obtain 



7(z) oc {z/R^Y'\ 



(34) 



This scaling was confirmed numerically as well [27, 

28]. 

It is necessary to stress that relation (32) takes 
place only if one can neglect the curvature of mag- 
netic surfaces. Indeed, for magnetically dominated 
case in the limit r J> rp the leading terms in 2D GS 
equation can be rewritten in the simple form [18] 

--n.V(B2)--^ + ^-— („.e.)=0. (35) 

Here i?c is the (poloidal) curvature radius of mag- 
netic surfaces, and n = V^'/|V4'|. Neglecting now 
the curvature term and using standard relations 
B^ « Bpr/iiL and B^ - w B^ 77.2^ resulting from 
(2) and (6), we return to (32). On the other hand, if 
the curvature is important, then one can neglect the 
first term in (35), and we obtain 

l^{Rjrf\ 



(36) 



This scaling taking place for split-monopole geom- 
etry outside the fast magnetosonic surface corre- 
sponds to 7 « 0-1/3 1^1/3 (^/^p) [29, 30]. Remem- 
ber that for r < rp we have "linear" acceleration 
(32). Thus, the effective particle acceleration can 
take place only if rjet ~ (tRi^, and if the curvature of 
magnetic surfaces is not important. 
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4.1.4. In the center of the self-similar domain 

The approach under consideration allows us match- 
ing the self-similar solution to the rotational axis. 
Indeed, for relativistic self-similar invariants 



(37) 
(38) 
(39) 
(40) 



the solution of the 2D GS equation for \1/ > \l/b has 
the form ^{p, 9) = p^/^@{6), where p is the spherical 
radius. Hence, far from the equatorial plane where 
z r {9 <^1, p K, z) one can write down 



*(p,^) =^pi/''r. 



(41) 



As a result, the cylindrical radius of the boundary 
\1/ = \l/b can be written as 



{Z) = ^-l/a,„l/a^l-l/afc_ 



(42) 



Let us consider now the central part of a flow 
^' < \E'b. If again 6* <C 1, one can integrate ID cylin- 
drical equations (6) and (10) considering 2; w p as 
a parameter. Assuming that Op = and ry = % 
for \1/ < \l/b and using solution (28) we have for 
Ml{z)=M\r^) 



Ml{z) = 



3-3/ab 



(43) 



On the other hand, for > \I'b (r > rb) one can seek 
the solution in a form A^^(r) = Mh^ir/rhY- As a 
result, equations (6), (10) give 



a = 2, e = 3 - 66. 
Substituting now r w z9, we obtain 

= C9\ 



(44) 



(45) 



The coefficient C oc M.h'{z)z^ /r^{z), in agreement 
with self-similar property, does not depend on z. 

4.2. Nonrelativistic flow 

4.2.1. Central core 

For nonrelativistic case in the central part of a flow 
one can use general expressions (9) 

(46) 



non-dimensional current io = j/jcj depending now 
on the angular velocity fip- For fip ^ f^cr, where 

1/2 

ficr = ^ ( 1 , (47) 



corresponding to particle dominated outflow near the 
star the 2D problem can be solved analytically [20, 

25], and we obtain io = c/fjn. For magnetically dom- 
inated flow near the origin one can write down [10] 

io « c/«i„ (^7FMr)"'/^ (48) 
Nevertheless, for < ^I'in, where 



io^l 



(49) 



the flow remains particle dominated: « uf„/2. 
Remember that in the nonrelativistic case the flow 
can pass smoothly the critical surfaces only if 
W^part(''F) ~ Wemirv)- Thus, the flow at large dis- 
tances is to be particle dominated. 

Solving now equations (22) and (23) for sub 
Alfvenic flow Mq < 1 we obtain that poloidal mag- 
netic field remains constant inside the jet up to the 
very boundary. Thus, one can put B{Q) — i?cxt- But 
such a fiow can exist only in the presence of large 
enough external magnetic field Bgxt > B{r-p)- For or- 
dinary YSOs B{rp) ^ 10~^ G, so sub Alfvenic flow 
in a jet cannot be realized. 

On the other hand, for super-Alfvenic cold out- 
flow one can find that the term oc in (23) plays no 
role. It means that here H « const, and we return 
to jet-like solution (1) [16, 20]. But, as was already 
stressed, in the presence of finite external pressure 
it is possible if the central core r < rcorc = ^in/^^ 
contains almost all magnetic flux \l/o. This can be re- 
alized only for slow rotation Qp <C 0,^. In this case 
magnetic field on the axis cannot be much smaller 
than Bmin = *o/7rr^ore: 

^' = HiT^J^y 

Accordingly, *core = *o/ln(l -|- Bmin/-Bext)- This 
structure was reproduced numerically as well [10]. 

But for fast rotation ^ ^cr the core mag- 
netic fiux \l/core is much smaller even than the flux 
^in (49) within the central part of a flow: 

^core , 



2cM2 



< 1. 



(51) 



It means that the cold cylindrical flow resulting from 
the interaction of fast rotating supersonic wind with 
the external media cannot be realized. 
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4.2.2. Heating in oblique shock 

To resolve this contradiction, one can propose that in 

real nonrelativistic jets an important role may play 
the finite temperature. E.g., the additional heating 
can be connected with the oblique shock near the 
base of a jet [31, 32]. It is well known that such a 
shock is needed to explain the emission lines observed 
in YSOs [33] .This situation is alike the hydrodynam- 
ical supersonic flow meeting the wall. This analogy 
is all the more reasonable as the non-relativistic su- 
personic outflow is to be particle dominated. 

To evaluate the thermal terms in equations (22) 
and (23) we consider pure hydrodynamical shock 
wave swifting spherically symmetric supersonic flow 
into cylindrical jet. Knowing the swifting angle, one 
can determine the entropy jump As as a function 
of particle flux, all other four invariants being the 
same as in front of a shock. As a result, we found 
that fastly rotating jet with Of ^ ^^cr heated in a 
shock is to have core-jet structure (30) with a < 2, 
(3 > 0. Hence, it can be realized in the presence of ex- 
ternal media. Obtained jet parameters (T ~ 10^ K, 
~ 10 km/s at r ~ 10 A.U.) [34] are in agreement 
with observational data. 

4.2.3. In the center of the self-similar domain 

The procedure similar to Sect. 4.1.4. for nonrelativis- 
tic particle dominated flow gives for self-similar re- 
gion (^' > *b, 6* < 1, oc *-''') that M"^ = Ce\ 
e = 2 — 46', as in numerical simulation [35]. 

5. Conclusion 

Thus, cylindrical GS equation has definite advan- 
tages in comparison with standard self-similar ones. 
Using this approach it was demonstrated that in rel- 
ativistic case effective particle acceleration can take 
place only if rjet ~ (tRi^, the curvature of magnetic 
surfaces playing no role. For nonrelativistic fiow we 
found that the heating in oblique shock near the 
base of a jet must play the leading role for mag- 
netically dominated flow. In both cases the magnetic 
flux within the central core was determined. 
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